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Abstract. Recently a new form of quantum memory has been proposed. The storage 
medium is an ensemble of electron spins, coupled to a stripline cavity and an ancillary 
readout system. Theoretical studies suggest that the system should be capable of storing 
numerous qubits within the ensemble, and an experimental proof-of -concept has already been 
performed. Here we show that this minimal architecture is not limited to storage but is in 
fact capable of full quantum processing by employing measurement-based entanglement. The 
technique appears to be remarkably robust against the anticipated dominant error types. The 
key enabling component, namely a readout technology that non-destructively determines "are 
there n photons in the cavity?", has already been realised experimentally. 
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1. Introduction 

Interactions concerning single spins are usually very weak, which imposes a huge challenge 
on the feasibility of their manipulation and measurement, and hence their potential for 
quantum information processing {QIP). Inspired by earlier work 03 13, it was suggested in 
2009 that electron spin ensembles on the other hand could exhibit strong couplings with a 
collective field via magnetic dipoles (see Fig. |l(b)D [31, exploiting the extremely small mode 
volumes in stripline cavities. Such strong couplings have also been recently demonstrated by 
independent experimental groups flUOl. 

Moreover, it was further proposed that the various spin waves in the ensemble could 
be used as independent quantum registers 0, for which a proof-of-principle experiment 
was carried out demonstrating the storage and readout of multiple classical excitations J6]|. 
Therefore, so-called holographic QIP schemes d [8l [9j should be possible by selectively 
coupling modes of the spin ensemble to the superconducting transmission line cavity J3]|. 
However, promoting a memory system to a full QIP device is non-trivial because the 
allowed operations in such a system are analogous to those in linear optics, and thus the 
provision of entangling gates requires non-linearity supplied by other components. One 
suggestion involves integrating a transmon Cooper pair box (CPB) ifTOl in the cavity; moving 
a memory qubit to the CPB would facilitate universal qubit operations, as well as the readout 
procedures 151171 [IT1. 

However, CPBs are not ideal for representing qubits due to their short coherence 
times IfTOl . Therefore in this paper we propose a composite system (see Fig. [[]) where such 
devices are only used for measurement, but in such a fashion that in addition to simple readout, 




(a) Circuit diagram of the integrated device for ensemble QIP. (b) Storage cavity coupled with an 

electron spin ensemble. 



Figure 1. (a) Circuit diagram of the integrated device for ensemble QIP, with figure adapted 
from fT3ll (see details of physical implementation on a chip in lil3l ). The number of photons 
in storage cavity affects the transmon state, which is indicated from microwave probing the 
measurement cavity; (b) Physical setup of the storage line cavity coupled with an ensemble of 
TV = 10 11 electron spins (N@Cgo) doped on substrate, with an average coupling g ~ 2ir x 20 
Hz. A bias field of 180 mT is required to bring the spin Larmor precession in resonance with 
the cavity at a resonance frequency of 2ir x 5 GHz (3J. A switchable linear magnetic gradient 
field for appropriate time lengths (gradient pulses) is required in order to access the different 
collective modes of the ensemble. Figure adapted from |3], with the CPB qubit removed. 
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also entangling projections can be implemented. We envisage the use of a system that is both 
number resolving lfT2ll and also non-destructive; recent experimental demonstrations of such 
systems have been accomplished lfT3l [T4| . The realisation consists of two cavities coupled 
via a transmon CPB; measurement of the number of single photons in the storage cavity is 
achieved via a read-out procedure involving the second cavity lfT3l . In effect this allows one to 
ask the question, are there exactly n photons in the storage cavity! If we wish to subsequently 
ask about a different n, we would employ a flux bias to tune the transmon frequency (a 
nanosecond timescale process). In the ideal case when no errors are present, if the result 
is NO then any coherent superposition of photon number states other than \n) is preserved in 
the storage cavity |[T3l . which is coupled to the electron spin ensemble 0. Appropriate bias 
field and magnetic gradient pulses are applied for resonantly accessing particular modes of 
the ensemble as discussed in Ref. 0. 

In contrast to previous proposals, here the CPB no longer performs the role of gate 
operations on the mode qubits nor to store the qubits. After a brief review of the basic 
physics in the spin ensemble we shall show that, by using dual-rail encoding, a universal 
set of quantum gates can be implemented for the logical qubits. Any single qubit rotation 
can be achieved by applying appropriate magnetic gradient pulses and adjusting the bias field 
when necessary. Importantly, the two-qubit parity projection that we shall presently discuss 
can enable general quantum computing through, for example, the creation of graph states |[T5i 

2. Modelling Collective Mode- Cavity Coupling 

Suppose the ensemble of N electron spins is in its ground state |0) = |0...0) and that the 
cavity contains a single microwave photon. If the bias field is such that the spins' Zeeman 
splitting is resonant with the microwave photon, then the ensemble will collectively absorb 
that photon on a timescale that is proportional to 1/ y/N. The ensemble state after the photon 
absorption is then given by 



where the sum is over all possible spin-flip ("1") positions q in the ensemble (see (31). Here, g q 
is the cavity coupling strength with the q th electron spin in the ensemble, and g is the average 
strength. The above state and state |0) together form an effectively closed two-level system 
(i.e., a mode qubit) Q. Now if no parameters are changed then the quantum of energy will 
"flip-flop" back and forth between the cavity and the collective state. 

However, if we wish to stop the flip-flopping, i.e., decouple the cavity from this memory 
mode, then we will temporarily apply a linear gradient in the magnetic field such that each 
spin will acquire phase at a different rate. If the gradient pulse causes the field to vary in the z 
direction, and writing the coordinate of the q th spin as z q , then the collective state becomes 




(1) 




(2) 
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Here, parameter k depends on the strength of the gradient and its duration £, and consequently 
the various terms in Eq.[2]have developed relative phases with one another 0. This prevents 
the y/N -enhanced mode-cavity coupling, and a single excitation is thus stored in the spin 
ensemble. Indeed, if we were to introduce a new photon into the cavity at this stage then that 
photon will resonantly transfer to the ensemble almost independently of the presence of the 
former excitation Q. After the application of the second gradient pulse, this procedure can 
be repeated yet again. 



As long as the number of the excitations n ^ N with an upper limit of n r 



< 



N 



(since the coupling enhancement of \J N — n ~ yN is still very large for small n), single 
excitations can be independently stored into and read out of the different collective modes 
i by appropriately applying impulses with ki = fc(£i) 0. The superradiant state |^i(0)) 
corresponds to the k = mode, and each time a fc;-piilse is applied it maps the k = mode to 
the i th mode. Therefore, to access a particular mode i a — fcj-pulse is applied to map it to the 
k = mode, which then interacts with the cavity strongly upon resonance. A fcj-pulse is then 
applied to map the mode back. To stop the the mode-cavity coupling, one can simply tune the 
bias field such that the cavity is out of resonance with the ensemble spins 0). 

Therefore, the Hamiltonian coupling a particular mode i with the cavity for a single 
excitation with energy e is 



= e (mjmj + c^c) + J(m\c + m^) 



e J 
J e 



(3) 



in the basis of |1m0c)> |0a/1c)- Here, m\ (c^) and m, (c) are the corresponding mode 
excitation (cavity photon) creation and annihilation operators, and J = y/Ng is the effective 
coupling strength being ~ 2n x 6 MHz for the parameters in Fig. |l(b)| This effective coupling 
J is orders of magnitude larger than the decay rates of collective spin excitations and the cavity 
decay rate BH. The corresponding time evolution operator for a single excitation is then 



-ict 



cos Jt 
-i sin Jt 



-i sin Jt 
cos Jt 



(4) 



After a time t = r := jj ~ 40 ns, a full SWAP of one single excitation has occurred between 
the cavity and mode i, whereas for t < r only a partial SWAP operation has taken place. 

When two or more excitations n exist between a particular mode % interacting with the 
cavity, the coupling strength and thus the time required for an exchange are adjusted by a 
factor of \fn. The cavity-mode Hamiltonian for two excitations is then 

/ 2e \/2J \ 



in the basis of |2m0c), 
S^\t)-- 



\ 



2e V2J 
V2J 2e V2J 
y/2J 2e 



(5) 



/ 



Wlc), |0m2c}> with the corresponding time evolution operator 

/ cos 2 Jt - -j= sin 2 Jt 
sin2Jt cos 2 Jt 

cos 2 Jt 



-2iet 



\ 



-sin 2 Jt — ^sin2Jt 



V2' 



sin 2 Jt \ 
/ 



-^sin2Jt 



(6) 
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Note that a full SWAP of two excitations still occurs after the same time t = t as for the single 
excitation case. 

For simplicity, from this point on whenever we say a particular mode i is interacting 
with the cavity for some time t, a — (+)/cj-pulse is by default applied right before (after) the 
interaction which has a duration of t with the bias field necessary for resonance. Note that the 
time for the necessary gradient pulses can be shortened by increasing the gradient strength 0. 



3. Single Qubit Operations in Dual-Rail Encoding 

We encode the qubits in the dual-rail representation, where each logical qubit occupies two 
collective modes of the spin ensemble (see Fig. El). A logical qubit \Q\)l = «|0)l + /3| l)x 
corresponds to the physical state of modes \MiM 2 )m = «|10)m + /3|01)a/, where the ket 
notation denotes the number of excitations in the relevant mode for the physical qubits. 
Adopting a dual-rail encoding of course doubles the resource cost. However, the capacity of 
the ensemble quantum memory is of the order of V~N ; assuming N = 10 11 , this is therefore 
not a practical constraint. 

An experiment would begin with cooling the spins so that there are a negligible number 
of excitations in the ensemble, which is then approximately in the ground state |0). For 
each logical qubit we then load the cavity with a single photon, and perform a SWAP gate 
between the cavity and the appropriate memory mode to represent |0)l- A logical qubit \Q) l 
of arbitrary state can then be prepared through an appropriate single qubit rotation on the 
Bloch sphere. 

Any single qubit rotation can be formed from a combination of rotations about two 
different axes lfT6ll . In Fig. [2] we show how rotations about the x-axis, and the z-axis can be 
performed in the dual-rail enconding. Both these rotations involve SWAP operations between 
the memory modes and an initially empty cavity. 
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(a) X-type rotations: Xgate when t = r 
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(b) Z rotations: phase gates 



Figure 2. Single qubit rotations in two independent directions, whose combinations constitute 
any single qubit operations for the dual-rail encoded qubits. Note that the cavity is assumed to 
be empty initially. 



3.1. X Rotations 

Suppose that our logical qubit is initially in the general state \Qi) — cos f \0} L + e 1 ^ sin || 1) L . 
We thus start with the initial state ( cos f 1 10) m + e 1 ^ sin ||01) A /) \0}c, where C denotes the 
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cavity state. After the first full SWAP, the physical state becomes 

cos^|00> M |l> c + e^sin^|01) M |0) c . (7) 

To implement a rotation around the x-axis on the Bloch sphere, a partial SWAP is now 
performed between the cavity and mode 2, resulting in the following state: 

cos^(cos^100) M |l) c -zsin^|01) M |0) c )+e^sin^(cos^|01) A/ |0) c -zsin^|00) M |l) c )(8) 
or equivalently, 

a|00) M |l) c + 6|01) A /|0) c (9) 

where 

9 9 
a = cos - cos 9' — ie 1 ^ sin - sin 9' 
2 2 

9 9 
b = —i cos - sin 9' + e 1 ^ sin - cos 9' 
2 2 

and 9' = Jt. When the second full SWAP has completed, the empty cavity decouples from the 
qubit \Qi) = a\0) L + b\l) L , ignoring a non-detectable global phase. An X gate is implemented 
in this way when t = r, i.e., when 9' = |. 



3.2. Z Rotations 

To obtain a phase gate, we begin in the same fashion and obtain the state in Eq. [7] However, 
we now exploit the fact that applying a different magnetic bias field does not affect the phase 
evolution of the component cos §|00)m|1)c (which still acquires phase at the same rate e), 
since the cavity photon energy is not affected by the magnetic field. On the other hand, a 
change 5B in the bias field shifts the spin ensemble mode energy to e + 5(5 B), which results 
in a phase evolution of e -1 ^" 4 "* 5 )' for the part sin | |01)m|0)c during the time t. Ignoring the 
undetectable global phase e~ let , after the second full SWAP the empty cavity again decouples 
from the qubit \Qi) = cos ||0)l + e l ^~ 5t ^ sin | \1) L . The relative phase accumulated e~ lH 
manifests itself as a phase gate for the logical qubit, where the phase can be controlled by the 
time t and the bias shift SB. 



4. The Measurement-assisted Entangling Scheme 

Having established how to perform single qubit operations, we now require an entangling 
operation in order to produce a universal set of qubit gates |[T6l . In this section, we show 
how certain types of measurement can facilitate such an entangling scheme (see Fig. |3]). In 
essence, we perform a measurement on the parity of two chosen logical qubits. If the parity is 
found to be odd ( |01)l or 1 10) i ), then we proceed with the protocol. Conversely, if the parity 
is found to be even ( 1 00) l or |11) l ), then we reject the state. This is therefore a probabilistic 
entangling operation, however, it is well established that such operations suffice for efficient 
universal quantum computing ifTTl . Our approach can be thought of as a filtering process; we 
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detect the state |00)l and reject it, and then similarly we detect and reject the state 1 1 1) L . Any 
state that passes through this process without rejection must be in the odd parity subspace 
{|01}l, 1 10) l}- Our protocol therefore consists of two basic blocks of operations, each ruling 
out the possibility of having two excitations in the relevant two modes while not destroying 
the quantum coherence between contributions to the state with different occupation numbers. 
As before, the cavity is emptied before the entangling gate operation, and it is tuned out of 
resonance with the ensemble modes during each measurement procedure. 
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Figure 3. The entangling operation for the dual-rail encoded qubits, implemented by a parity 
projection. A 'building block' consists of all the operations inside the dashed box, and the 
parity projection requires two of such blocks. Each building block rules out the possibility 
of two excitations in the relevant two modes while preserving coherence between states with 
different occupation numbers. This requires two NO results for each block to the question 
whether there are 2 photons in the cavity or not. The initial cavity state is assumed to be 



empty, and the interaction times are chosen as follows: ii = jj 

and t 3 = 2r — ti — t 2 



20 ns, t 2 = ~ 40 ns 



fj ~ 20 ns. 



Let us now consider the first block of operations, concerning modes Mi and M 3 with the 
cavity. For each correctly encoded logical qubit, the two dual-rail modes together only contain 
a single excitation. However, for the M 3 mode-cavity interaction there is still the possibility 
of exchange of two excitations, originally from different qubits. This two-excitation exchange 
is governed by the Hamiltonian in Eq.[5l which has three distinct eigenvalues 2e, 2(e — J) 
and 2(e + J), with corresponding eigenvectors 

\vo) MO = -^(|20) - |02» 
|t;_) A/c = ^(|20)-v / 2|ll) + |02)) 

I^+)a/c = ^(|20) + v / 2|11) + |02)) (10) 

respectively. Therefore, if the M 3 mode-cavity state after the first full SWAP is |11) = 
-j^(\ v +) — \ v -))> ^ evolves to -L(|t> + ) + = ^=(|20) + 1 02) ) for an interaction time 

ti = fj ~ 20 ns. Performing a measurement designed to detect whether there are 2 photons in 
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the cavity |fT3l then removes the 1 02) component if the measurement result is a NO outcome. 
After a further exchange of t 2 = fj — 40 ns, the part 1 20) evolves to 1 02) which can be 
ruled out by a second measurement if the outcome is NO again. Thus, up to now when 
successful, we are assured that there will be at most one excitation between the modes Mi 
and M 3 , after it has been swapped back from the cavity. So far, we have only considered 
the two-excitation exchange, while there is also the exchange of one excitation. In order for 
the block to perform an identity operation, i.e., two SWAPs, on the one-excitation exchange 
subspace, a further interaction between the mode M 3 and cavity is required for a time of 
t 3 = 2t — fx — t 2 = jj ~ 20 ns, followed by a full SWAP between the mode Mi and cavity. 

Due to symmetry, the second block of operations has the same effect on modes M 2 and 
M4, ensuring that there is also at most one excitation in between, if successful. Hence, starting 
with two correctly encoded qubits \QiQ 2 ) L = («i|0)l + /3i|l>x) 0(a 2 |O)i + f3 2 \l) L ) or 
physically |M 1 M 2 M 3 M 4 ) = a 1 a 2 |1010) + ai/3 2 |1001) + /3i« 2 |0110)+/3i^2|0101), successful 
implementation of the entangling gate leaves us with the mode state A/"(ai/3 2 |1001) + 
/3ia 2 |0110)) or A/"(ai/3 2 |0 Ql lQ 2 ) L + /3ia 2 |l Ql 0Q 2 ) L ) in the logical basis, where J\f is a 
normalisation constant. This is an odd-parity projection operation, which generates a 
maximally entangled Bell state when we start with \+)q 1 \+)q 2 as the initial state of the logical 
qubits. Notably, such parity projections together with single qubit operations are adequate for 
universal QIP IfTTI 

5. Errors 

With the above set of universal quantum gates, our integrated device is capable of performing 
probabilistic QIP within the electron spin ensemble, in the ideal case when no errors are 
present. We shall now consider the various error sources, namely, collective excitation decay, 
imprecise timing in cavity-mode swap, cavity photon leakage, and measurement errors. 

As mentioned before, the decay rates of the collective spin excitations and that of 
cavity photon leakage are orders of magnitude smaller than the collective coupling J 
Moreover, errors resulted from imperfect pulse timings are all linear, and timing precision 
with a nanosecond resolution is feasible Q. Therefore, these errors shall be treated negligible 
as compared to measurement errors lfT3ll . which shall be evaluated for the entangling gate in 
this section. We start with the logical qubit state 

\+)qA+)q 2 = |(|ioio> + liooi) + 10110) + |0101» A , (11) 

with the cavity being \C) = |0). There are two types of measurement errors to consider, 
corresponding to different incorrect outcomes to the question whether there are two photons 
in the cavity or not lfi~3ll . 

5.1. Type I Measurement Error 

This type of error occurs when the measurement outcome gives a NO while the cavity is 
populated with two photons, and this occurs with an independent probability of 771 for each 
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measurement. It has an effect on the overall fidelity F of this entangling gate, reducing the 
degree of entanglement generated in the final state. With the cavity state |0)c decoupled from 
the modes at the end of the protocol, the final density state for the modes IM1M2M3M4) in 
this case is (see Appendix) 



p = _^_( |flr+)(tf+| + + a 2 ) ) (12) 



where 



i)iV " 1 2 



|*+) = ^=(|1001> + 10110)) — ► -^(|01> + |10)) QlQ2 (13) 
^ = 1 ((10020) -I2000))((0020|- (20001) + |ioio)(io1(j| 

a 2 = 1 ((10002) -I0200))((0002|- (02001) + ^ ^ 

In the ideal case where r\\ = 0, the resulting state after the entangling operation is 
the maximally entangled state \^ + ), as expected. For finite rjx errors, the final mixed state 
p describes that for bipartite quinits Q\ and Q 2 , where each quinit is a five-level system. 
Mapping the different physical states of the dual-rail modes |10), |01), |20), 1 02) and 1 00) to 
the logical quinit levels |0), |1), |2), |3) and |4), respectively for each Q { , Eq. [T4lbecomes 

gl ^l(" 42 >-' 24 » 2 '< 42 '-< 24 " + |00)(00| 

2\ 2 / Q1Q2 

Although also being entangled in the levels |2), |3) and |4), we shall only be interested in the 
entanglement between the two levels |0) and |1), which span the computational basis space 
for conventional QIP tasks. Fig. |4(a)| plots the fidelity F |[T8l of the resulting state p (as given 
by Eq.[T2l with the ideal state |\E' + ) (\E' + | against the error 771, 



F(p, = Tr (yfy/p\* + )(9 + \y/p) = ^jy^ (16) 

where |^ + ) is embedded in the five-level system by filling the levels that lie outside the 
computational basis with zero population. 

We see from Fig. |4(a)| that a high fidelity is possible as long as the error 7]i remains small, 
however, this does not automatically guarantee the presence of entanglement for the bipartite 
quinits. We therefore apply an entanglement measure to see how the degree of entanglement 
changes as rji increases. The simplest (but questionable) way of doing so is to apply Wootter's 
entanglement of formation E F |fi~9l to the state p projected onto the computational subspace, 
while ignoring the population in the other levels (see Fig. |4(b)D . The reason we consider this 
metric incorrect is that artificial entanglement is created by this projection process. Moreover, 
even small amounts of population in levels outside the computational subspace may result in a 
significant reduction in the degree of entanglement for the larger system in certain scenarios, 
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Figure 4. Plots for the final state p after the entangling operation, against the type I error r)i for 
each measurement, showing (a) the fidelity F with respect to \ ^ + ) (^ + |; (b) the Entanglement 
of Formation Ep when the population outside the computational subspace are simply ignored 
and projected out. 



as has recently been explored by Ref. Il20l . Thus instead of the entanglement of formation, a 
more useful measure involving the whole five levels of the bipartite quinits is needed. 

We apply the generalised m-concurrence measure in d dimensions as defined in [|2T| 
and derived in [|22l . where m denotes the number of parties involved and equals 2 in this 



case. In d dimensions, Cf ranges from zero for separable states to 2 ^~ 1 - > for the maximally 
entangled generalised Bell states, such as ^= J2i=o N*)- Moreover, it reduces to the Wootter's 
concurrence squared C 2 |[T9ll when d — 2. Fig. |5(a)| illustrates the 2-concurrence {dashed 
green) of the final state p (as in Eq. [121 with d = 5 against the error rji. This is marginally 
lower than the projectively implemented Wootter's concurrence squared (not shown here), 
for each finite rji ^ 0. The difference is insignificant, partly due to the fact that the state p 
in Eq. [12] is also entangled outside the computational space through the error components 
o\ and o"2 (Eq. [T51) . Interestingly, even when 771 = 1, the 2-concurrence is still greater 
than zero, although no information is learnt at any of the four measurements. This can 
be understood as follows: we have assumed that the phase between the measured state 
and the rest of the superposition is lost through the interaction between the system and 
the measurement apparatus (see Appendix) |[T3l . This corruption of phase information 
affects the subsequent evolution of the state and effectively causes a leakage of population 
from the components |00)q 1 q 2 (|H)q 1 q 2 ) in the computational subspace to |42(3))q 1 q 2 and 
|2(3)4)q 1 q 2 . Entanglement then appears to exist - and is registered by the 2-concurrence - 
both within and outside the computational subspace. However, we stress that this apparent 
entanglement is useless for our purpose due to a complete lack of information about heralded 
failures, which prevents the required post-selection of successful runs of the protocol. 

As discussed above, our final state generally contains both useful and also 'spurious' 
entanglement. Therefore, we require a new metric for filtering out the latter kind to obtain 
a lower bound on the degree of useful entanglement. We shall achieve this by identifying a 
suitable theoretical operation that can be applied to the final state. It turns out to be insufficient 
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Figure 5. (a) Plots of the generalised 2-concurrence Cj in 5 dimensions against the r}\ error, 
for p (dashed green), p with dephased population outside the computational subspace {solid 
blue), and p with all population discounted into the computational subspace (solid red); (b) 
Entanglement of formation Ep plot for the discounted final state pd (red in (a)) against rji. 



to replace all population outside the computational basis with a completely mixed state in 
these levels, since the resulting curve (solid blue in Fig. |5(a)D still does not hit zero as 771 — > 1. 
(In this case only the entanglement outside the computational subspace has been removed, 
but the entanglement inside of it is still present. Since no information has been learnt in the 
measurement, this entanglement is of the 'spurious' kind which we want to discount.) A better 
approach therefore proceeds as follows: we add any population that has leaked from the levels 
\00)q 1 q 2 or I 11)q 1 q 2 back into computational subspace in a way that mixes the resulting state 
to the largest possible extent: 

p d = -^—(\* + }(* + \ + f (|00)(00| + |11)(11|)) (17) 
I + 771V 2 /Q1Q2 

The Cf curve in 5 dimensions for this state then coincides with the square of Wootter's 
concurrence C 2 for two 2-dimensional systems, hitting zero as r/i — > 1 (solid red in Fig. |5(a)D - 
Indeed, this is what we would reasonably expect for maximal error rji when the measurement 
fails to report any useful information. Since only two levels are involved for p d , we can also 
plot its entanglement of formation against 771 in Fig. |5(b)[ which is lower than that shown in 
Fig. |4(b)| as expected. We see that in order to achieve a high degree of useful entanglement as 
the result of the parity projection operation, 7/1 has to be small. For example, Ep(pd) — 0.75 
when 771 = 10% which is achievable with the current technology lfT3l . 

5.2. Type II Measure Error 

This type of error occurs when a measurement outcome reports a YES while there are not 
two photons in the cavity, and this happens with an independent probability of r/ 2 for each 
measurement. For this kind of error, the entangling procedure is STOPPED and we have to 
start all over since a YES is reported. Consequently, this only affects the overall probability 
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that we reach the end of the protocol without a heralded failure. We found this by considering 
the overall probability of passing all four measurement tests, and thus obtain 

D _ 7(1 - ?7 2 ) + 771 (6 + T7i)(l- 77 2 )+ 77i (5 + 2r 7l )(l - 773) +77! (4 + 3770(1 -772) +77! 



; " 8 7 + r/i 

for the entire entangling operation (see Fig.©. 



Ps 




(a) 771 = and 0.1 



6 + 2t/! 5 + 3771 



Ps 




(b) 772 = 



Figure 6. The probability Pf (a) for a type II error 772 with 7/1=0 (blue) and 771 =0.1 (red), 
respectively; (b) for a type I error 771 with 772 = 0. In the latter case, note that in spite of the 
increase in the probability Pf the final state is still degraded as 771 becomes larger. 

In the absence of any error, the probability Pf of reaching the end of the protocol is 50%. 
Pf decreases with r/ 2 for fixed 771, and increases with 7]i for fixed 772, although the quality of 
entanglement is of course reduced in this case. For the 10% error in 771 which gives rise to an 
80% concurrence, a type II error of 772 = 10% lfi~3l results in a probability Pf ~ 36.4%, i.e., 
roughly a third of the times. 

5.3. Cavity Photon Leakage 

Measurement errors are only important during the entangling scheme, while all other errors 
also affect the single qubit gate fidelities. Cavity photon loss is likely the next dominant error 
source, whereas timing errors and collective excitation decay can be considered negligible in 
comparison. 

For this reason, we shall now briefly discuss an adaptation of our proposed circuit which 
can achieve a reduction of the cavity leakage effect. For each measurement procedure a 
Control-NOT (Cnot) gate is first implemented on the transmon, which takes about 50 ns. 
The Cnot gate is controlled by the number of photons in the cavity, so in this case the Cnot 
gate is triggered whenever there are two photons in the cavity. The measurement cavity is 
then read out using a microwave probe, taking a time of about 400 ns lfl~3l . 

As the probe measurement thus takes a much longer time than the Cnot gates and the 
SWAP operations, we can adapt our earlier circuit to reduce the number of read-out steps 
required in each block of the entangling operation (see Fig. [7]). This modified block achieves 
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Figure 7. The modified circuit for the whole integrated device during the entangling operation, 
where the number of measurements is halved. Nonetheless, the circuit still consists the same 
blocks of basic operations. 



the same effect as our previous circuit, i.e. it removes all components of the wavefunction with 
double excitations in the relevant two modes while not disturbing the remaining components 
of the state. For a cavity with photon decay time of P ~ 3.2/is |[T3ll we expect a photon 
leakage of about 15.6% during a 0.5 /xs operation, but this source of error can be further 
suppressed by increasing the Q-factor of the storage cavity in the integrated device. 



6. Readout Scheme 



Full quantum state tomography can be carried out on the logical qubits to determine the 
dual-rail mode states. This is achieved through a sequence of SWAP operations followed 
by measurements on the storage cavity. This can also be extended to consider the entangled 
bipartite quinits as in Eq.[[2l by measuring the number of photons in the cavity swapped from 
each mode: Three measurements are then needed for each mode to determine whether there 
are zero, one, or two photons in the cavity. 

While the above-described full state tomography is straightforward to implement in 
principle, the number of measurements required for each run of the circuit might be 
challenging in a first experimental demonstration. Therefore, we propose the following 
alternative as a simpler verification of the entangling protocol: the aim is to use appropriate 
measurement statistics on the final state to reconstruct a density matrix that is consistent with 
the one given by Eq. [FT] To achieve this, we start with repeated runs of the solid circuit in 
Fig. [Ma), with each measurement only checking whether there is a single photon in the cavity 
or not. An answer of NO means that there must have been either zero or two photons in the 
relevant mode Mj. Our default assumption in this case will be that the mode was in the state 
\0}m v Only when we obtain the total combination 1 0000) ^a^MsM^ do we assume that there 
must have existed two excitations in one of the four modes. However, such states only occur 
due to errors in the protocol; they lie outside the computational basis and do not contribute 
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towards the final entangled resource state. As discussed in Sec. |5.1| we account for these errors 
by adding their population to the computational subspace in a way that ensures that the quality 
of the desired resource is not artificially boosted. Performing many runs of this circuit then 
allows us to fill the diagonal entries of the density matrix (in the {|0)Q i , 11)^} basis), and to 
compare them with Eq. \T7\ 



Ml (Qla) ■ 

M2 (Qib). 

Ms (Q2a) ■ 
M4 (Q2b) ■ 



Cavity 



-4- 



S(T) 



Measl 



Meas2 



) [SCto)l 



4- 



S(T) 



Meas3 



Meas4 



Figure 8. Qubit measurement in two independent bases: (a) with only the solid ciruit, 
corresponding to measurement in the logical basis of { 1 0} , |1)q;}; (b) with the full circuit 
including the dashed gates, corresponding to measurement in the logical basis of 
|-)qJ for t = |. 



However, the above procedure alone is not sufficient to distinguish quantum 
entanglement from classical correlations. We therefore also carry out repeated runs of the 
complete circuit including the dashed gates in Fig. [8Kb), corresponding to measurement in the 
logical basis of {|+)g t , |— )q 1 } for t — |. This enables us to obtain the coherences of the 
density matrix, which can then once more be compared with Eq. [FT] Note that in measuring 
the coherences, population outside the computational basis is no longer relevant, and can be 
directly discarded. Combining the two similar procedures described above, one can verify 
whether the entangled resource at the end of the protocol is consistent with the discounted 
qubit state pd for a known fixed rji error, up to a certain confidence level. 

7. Conclusion 

We have devised measurement-based entanglement protocol for dual-rail encoded qubits in 
the collective excitations of an electron spin ensemble, coupled to a superconducting resonator 
linked with a measurement apparatus. An detailed analysis of the predominant error sources 
of the integrated device indicates that probabilistic QIP with the spin ensemble should be 
feasible with the current technology. 
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Appendix 

In this Appendix, we go through the evolution of the total mode-cavity state under the 
entangling operation with type I measurement error r q 1 , to derive the final state p as in Eq.fTH 
Since there are two excitations in the whole system, the phase evolution e ~ 2let is global and 
ignored at all times. 

The initial total state \M\M 2 M- i M^C) of the modes (Eq. fTT]) and the cavity evolves to 

if 100200) + 100002) , |oooii) -jjooiio) | |oiioo) - ijoiooi) | |0101Q) \ (18) 



2V y/2 y/2 V2 

after the first full AfiC SWAP followed by the M 3 C coupling S(tx). When the first 
measurement gives an outcome of NO, we are left with the density state 



Pi 



( y -\A 1 )(A 1 \ + -77i|00002>(00002|) (19) 



7 + ?7i 
where 

]Al) = ./|(^|oo 2 qo )+ -'I° w11 > - B + B -;l 01 ° 01 ) + |o 101 o)) ( 20) 

V 7 V y/2 \/2 v 2 ' 

all normalised. Note that here with the set-up in lfT3l . we have assumed that the measurement 
apparatus interacts with the cavity-ensemble system regardless of its ability in reporting 
correct answers. In so doing, even when r\\ — 1 where no information is to be learned from 
the measurement result, the state has still changed, namely, dephased according to Eq. [T9l 

The M 3 C coupling S(t 2 ) is a full SWAP for both the exchange of one and two excitations. 
Conditional on the second measurement giving an answer NO as well, the density state 
becomes 

P2 = 7^ A 6 \ A 2)(M +r/! 1 00002) (00002 1 + 1 00200) (00200 1 ) (21) 

2(3 + 77!) \ / 

where 

\A 2 ) = — ((—(00110) +i|00011» + (— i|01001) - |01100)) + >/2|01010>) (22) 
all normalised. 

The remaining operations in the first block results in the density state 

^\A h )(A fl \ +71lPo ) (23) 



where 



and 



3 + 77! 

\A h ) = -^=(J10010) - |01100) + |01010)), (24) 

1 / (100200) - |20000))((00200| - (200001) A 
Po=2{ 2 - + |10100)(10100|J (25) 

which does not evolve further during the rest of this entangling gate operation since the second 
block of operations never access the modes Mi and M 3 . At this point, the cavity is empty and 
decoupled from the modes. 
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Similarly due to symmetry, upon the completion of the second block of operations on the 
modes M 2 and M 4 , the cavity state \Q)c is again decoupled from the modes and we have for 
the modes M X M 2 M Z M^ the density state (Eg. [HlfT?]) 



where 



l*+) = -L(- |iooi> -|oiio)) 

G\ = 1 ((10020) -|2000)K(0020|- (20001) + | 1010)(1010| ) 
ff2 = 1 ((10002) -|0200)K(0002|- (02001) + |Mfll)<0101| ) 
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